Solitons in the Calogero model for distinguishable particles 
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Abstract 

We consider a large — N, two-family Calogero model in the Hamiltonian, collective- 
field approach. The Bogomol'nyi limit appears and the corresponding solutions are 
given by the static-soliton configurations. Solitons from different families are local- 
ized at the same place. They behave like a paired hole and lump on the top of the 
uniform vacuum condensates, depending on the values of the coupling strengths. 
When the number of particles in the first family is much larger than that of the 
second family, the hole solution goes to the vortex profile already found in the one- 
family Calogero model. 
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The Calogero system is a class of exactly solvable models in one dimension [1-3]. 
The model has found wide applications in areas as diverse as condensed matter [4], 
black- hole physics [5] and two-dimensional string theory [6]. The ordinary Calogero 
model describes N indistinguishable particles on the line which interact through 
an inverse-square two-body interaction. As far as the distinguishable particles are 
concerned there are a few generalizations of the Calogero model to models of particles 
with different masses and with -j- couplings depending on the labelling of the 
particles coupled [7-14]. A multispecies one- dimensional Calogero model with two- 
and three-body interactions was treated in the SU (1,1) algebraic approach in 
[12, 13], while its matrix formulation was presented in [14]. Although it was possible 
to find an infinite number of exact eigenstates and eigenenergies, the set is not 
complete. 

Recently, it was shown in Ref. [15] that a natural supersymmetric extension 
of the Calogero model resulted in two-family Calogero models which are exactly 
solvable in some special cases. 

In an attempt to better understand the nature of the multispecies Calogero 
model, in this paper we transform the two-family model to collective fields. This 
transformation gives much more insight into the non-perturbative, solitonic sector 
of the theory. 

Let us start with the Calogero Hamiltonian describing two different families of 
particles in interaction [13]: 
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The first family contains N\ particles with mass mi at positions Xi, i = 
l,...,Ni, while the second one contains N 2 particles with mass m 2 at positions 
x a , a = 1,...,N 2 . The particles of the same kind interact and the corresponding 
coupling constants within each family are given by A and v, respectively. The 
particles of different kind also interact and the interaction strength between the first 
and the second family is denoted by k. We consider the parameters A, v and k 
positive. 

The Hamiltonian (1) describes the simplest multispecies Calogero model for par- 
ticles on the line, interacting with the two- and three-body potentials. Setting 
A = v = k and m\ = m 2 , we recover the ordinary iV— body Calogero model. 
The three-body terms in (1) trivially vanish in this case. General conditions for the 
absence of three-body interactions are given in [12, 13]. In the following we do not 
use any confining potentials. 

We can perform the similarity transformation 

H - U-'U-'U^HU^Un (2) 
to obtain a simpler but non-hermitian Hamiltonian 
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where the two- and three-body interactions have simply disappeared. The Jastrow 
prefactors are given by 

n A = n ( x i ~ x i) A ' 

i<j 
N 2 

a</3 

n K = n - X *) K (4) 

i,a 

and incorporate the conditions that the wave functions go to zero whenever the 
particles approach each other. The usual approach to the quantum-mechanical 
problem is to solve the eigenvalue problem with the Hamiltonian (1) using symmetric 
and antisymmetric wave functions, depending on the underlying statistics of the 
identical particles. Instead, we develop a collective-field theory of this system in the 
large — Ni and large — N 2 sectors of the Hilbert space. 

The collective-field theory for the two-family Calogero model is obtained by 
changing variables from the particle coordinates Xi and x a to the density fields 
p{x) and p{x) defined as 

Ni 

p(x) =J2$( X ~ x i), (5) 
i=i 

N 2 

p(x) = ^5(x - x a ). (6) 

a=l 

Such a change of variables is meaningful only if the particle numbers JVi and N 2 
go to infinity [16, 17, 18]. The Hamiltonian (3) can be expressed entirely in terms 
of p(x), p(x) and their canonical conjugates 



satisfying the following equal-time commutation relations: 

\p(x),ir(y)]=i5(x-y), (9) 

\p(x),n(y)]=i6(x-y), (10) 
\p(x),p(y)] = [ir(x),t(y)] = 0. (11) 
After the change of the variables, the Hamiltonian (3) takes the form 
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where j- denotes Cauchy's principal value of the integral. This Hamiltonian is 
still non-hermitian owing to the imaginary terms. In order to obtain the hermitian 
Hamiltonian, we have to rescale Schrodinger's wave functions of the original Hamil- 
tonian by using the Jacobian of the transformation from {x iy x a } to {p(x), p(x)}, 
as was suggested in Ref. [16]. After performing a straightforward algebra, we find 
the Jacobian J 

InJ = (1 — A) J dxp(x)lnp(x) + (1 — u) J dxp(x)lnp(x) 

—A J dxdyp{x)ln\x — y\p(y) — v j dxdyp{x)ln\x — y\p{y) 

—2k J dxdyp{x)ln\x — y\p{y). (13) 
The hermitian Hamiltonian is finally given by 
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where P stands for the principal part. The two terms, quadratic in the conjugate 

momentum operators 7r and fc, represent the kinetic energy of the system. The 

rest emerges as a quantum collective-field potential. The last terms are singular and 

do not give a contribution in the leading order in the and expansions. 

They should be cancelled by the infinite zero-point energy of the collective fields p 

and p. 

To find the ground-state energy of our system, we assume that the corresponding 
densities are static. Since their momenta are vanishing, the leading part of the 
Hamiltonian in the and expansions is given by the effective potential 
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Its form makes the Bogomol'nyi bound apparent. The potential is positive semi- 
definite and its contribution to the ground-state energy vanishes if there exist posi- 
tive solutions of the coupled equations 

A ~ 1 d x p + A r dypjy) + f dypjy) = Q 
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It is evident that there always exist uniform solutions 

p(x)=p , p(x) = p . (17) 

We have not been able to obtain analytic solutions to these equations for any 
values of the parameters A, v and k. However, if we further simplify our model by 
the assumption that there are no three-body interactions between the particles in 
the starting Hamiltonian (1), we obtain the conditions ( see Ref. [13] ) 



k 2 = Xu, 
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In this particular case, the solutions of equations (16) are always interrelated by 

p vh ' v ~ h ~ p xh ~ x ~ h . (19) 

Note that for A = u, the condition (19) implies proportionality between p and p 
and this means that we are dealing with the one-family model. 

If we further assume that k — 1 ( weak-strong coupling duality in Ref. [19] ), 
we end up with the condition 

pp = c, (20) 

where c is some positive constant. This very condition allows us to find new soliton 
solutions to the coupled equations (16). In fact, there is only one relevant equation, 
let us say for p, for example 

x ~ ldxP | \ f dypjy) { f dy _ Q 

2 p 7 x - y 7 p(y)(x - y) 



This equation can be solved by a rational Ansatz 



x z + tr 



where a and 6 are positive constants. By using the Hilbert transform 

dy 1 7r £ 



we find the conditions 



x — yy 2 + a 2 a x 2 + a 2 ' 
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(23) 



l-A + ^V-5 2 ) = 0. (24) 



The soliton solution for the first family is given by 



P(X)= * + Vs'V (25) 
while the solution for the second family looks like 

~t \ c c 1 — A a ,^„ s 
p(x) po n x 2 + a 2 

We note that both solutions are localized at x — 0. For large values of x, the 
soliton solutions approach the uniform solutions found before. This yields one more 
condition 

P(oo) = Ni = P^_ (2?) 

p(oo) N 2 c ' 1 J 

which, together with the conditions (24), finally fixes the parameters a,b and c, 
namely 

JVx(A-l) 

a - 
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(1-A) 
Ap 7r(l- A^/A 2 ^ 2 )' 

c=*'* (28) 
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The particle number of the first soliton is 

A - 1 



J dx(p(x)-po) = ^-, (29) 

while that of the second soliton is 

J dx{p{x) -po) = l-X, (30) 

where po — It is interesting to observe that these numbers are generally not 
integers. 

For A < 1, the first soliton behaves like the hole in the condensate p and 
the second one behaves like the particle above the condensate pQ. The roles are 
interchanged for A > 1. It is interesting to study the limiting case c — > 0. This 
can be achieved only if Ni is much larger than N 2 . The condition (24) gives 
a — > (the case b — > oo destroys the solution). The first-family soliton reduces to 
the "vortex" profile 

* x) = Po 1?Tp> b = ^ (31) 

The above solution is called a vortex because the density goes to zero at x — 0. It 
has already been obtained in the collective-field approach to the one-family Calogero 
model [18, 20]. The second-family lump solution transforms into the sharp delta- 
function profile 

p(x) = (1 - X)5(x). (32) 
In deriving this result we have used the Lorentzian representation of a delta function 

6(x) = lim --=-^-2. (33) 

We have thus demonstrated that our solitons have vanishing energy in the leading 
approximation. It does not cost any energy to create such a pair of correlated 
solitons. 
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What about the other solutions to the coupled equations (16)? The authors of 
Ref. [19] have claimed that there exists a multi-vortex solution in the first family 
accompanied by the sum of the delta-function profiles for the second family. The 
zeros of the p(x) (i.e. the positions of the vortices) are simultaneously the points 
at which p(x) diverges (i.e. the positions of the particles). According to Ref. 
[19], "the collective field p(x) describes solitons as particles, whereas p(x) gives 
a microscopic description in terms of fields". This is the essence of the so-called 
particle-vortex duality. We do not agree with this interpretation since p(x) and 
p(x) describe two different physical systems and, consequently, cannot be applied 
simultaneously to the same object ( soliton in this particular case). Moreover, we 
doubt that there exist multi-vortex solutions at least not of the form suggested in 
Ref. [19]! Nevertheless, let us investigate this possibility more carefully in the case 
of the hypothetical two- vortex solution. This solution can be reached by the limiting 
procedure (c — > 0), starting from the two-hole solution in the condensate p 



where a, a and b, b are complex parameters. This Ansatz for p(x) takes 
into account the positiveness of the particle density automatically. Furthermore, 
the rational function of two polynomials with the fourth degree in x in principle 
ensures the two-hole ( and for P—- p i the two-lump) structure. 

Inserting the expression (34) for p and - for p in equations (16), and by 
using the Hilbert transform 



where the + sign ( — sign) is taken for a in the upper (lower) half of the complex 



p(x) = po 



(x 2 — a 2 )(x 2 — a 2 ) 
{x 2 -b 2 ){x 2 -6 2 )' 



(34) 




(35) 
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plane, we obtain the following system of algebraic equations: 

A - 1 ± (a 2 - b 2 ) — — — = 0, 

Poa cr — cr 



1 - A ± —^(6 - a = 0. (36) 



Let us set 



a = a + i/3, b = 7 + i5, (37) 

where 7 and 5 are real numbers. The real and the imaginary parts of the 

complex equations (36) are, respectively, 

A-l±^(a 2 -/? 2 - 7 2 + 5 2 ) = 0, 
PoP 

1 _ X ±^( 1 2 -5 2 -a 2 + f3 2 ) = 0, 
Wf3 2 - 4 7 2 <5 2 + (a 2 - P 2 - 7 2 + 5 2 )(3a 2 + /3 2 + 7 2 - 5 2 ) = 0, 
4 7 V - Aa 2 (3 2 + ( 7 2 - 5 2 - a 2 + /3 2 )(3 7 2 + 6 2 + a 2 - (3 2 ) = 0. (38) 
The only solutions are 

a = 7 = 0, /5 2 = <5 2 , 

a 2 = 7 2 , = 5 = 0, 

a = / 3 = 7 = 5 = 0. (39) 

This brings us back to the uniform solution p(x) = p . Similarly, one can show that 
the Ansatz p(x) = Po ^ 2 ^j^2^^) ; for real parameters a\,bi,(i2 and 6 2 , leads to 
the one hole- lump pair solution, Eqs. (25) and (26). 

In conclusion, the only duality that can be considered is the trivial self-duality 
of the collective Hamiltonian (14). Obviously, it is invariant under the interchange 
of the two families, namely p <-> p, n <-> n, A <-> v and m\ <-> m 2 . 
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Our final remark is that our two-family Calogero model without the three-body 
interaction (18) can be viewed as the one-family model, but only in the leading 
approximation. Namely, comparing the effective potential (15) with the effective 
potential for the one-family Calogero model, one finds that they can be identified. 
To this end, we can define the new, effective one-family collective field p e ff = 
( ^" mi ) 1 ^ 3 (Ap + up), describing one-family particles with the effective mass m and 
effective one-family coupling strength p. Note that this identification cannot fix the 
values of the afore-mentioned parameters. However, owing to the presence of the 
d x p/p and d x p/p terms, which are suppressed, respectively, by the factors and 

with respect to the leading terms, this equivalence breaks down. The derivative 
terms are crucial for the appearance of the soliton solutions. In this respect, the two- 
family models display the new, paired solitons which do not exist in the one-family 
Calogero models. 

Although we have here focused on a two-family Calogero model, the discussion 
in this paper should be relevant to two- matrix models, since it is known that the 
Calogero model actually corresponds to the O(N), U(N) and Sp(N) invariant 
matrix models for A = §,1,2, respectively [21]. Our results can also be eas- 
ily extended to the models with more than two distinct families of distinguishable 
particles. The open problems that still remain are the question of the quantum sta- 
bility of the semi-classical solutions (25) and (26), and the existence of the possible 
moving-soliton solutions. We hope to study these issues in the near future. 

Acknowledgment 

This work was supported by the Ministry of Science and Technology of the Republic 
of Croatia under contract No. 0098003. 

12 



References 

[1] F.Calogero, J.Math.Phys. 10 (1969) 2191, 2197; 
F.Calogero, J.Math.Phys. 12 (1971) 419. 

[2] M. A. Olshanetsky and A. M. Perelomov, ITEP-90-1979 

[3] B. Sutherland, Phys. Rev. A 5 (1972) 1372. 

[4] F. D. M. Haldane, Phys. Rev. Lett. 60 (1988) 635; 
B.B.Shastry, Phys. Rev. Lett. 60 (1988) 639; 
H. Azuma, S. Iso, Phys. Lett. B 331 (1994) 107. 

[5] G. W. Gibbons and P. K. Townsend, Phys. Lett. B 454 (1999) 187 [arXiv:hep- 
th/9812034]; 

D.Birmingham, K.S.Gupta, S.Sen, Phys. Lett. B 505 (2001) 191. 

[6] J. Gomis and A. Kapustin, JHEP 0406 (2004) 002 [arXiv:hep-th/0310195]. 

[7] P.Forrester, J. Phys. A: Math. Gen. 25 (1992) L607. 

[8] C. Furtlehner and S. Ouvry, arXiv:hep-th/9407107. 

[9] D. Sen, Nucl. Phys. B 479 (1996) 554 [arXiv:cond-mat/9512014]. 
[10] A. Dasnieres de Veigy, arXiv:hep-th/9603050. 
[11] S.Maskhevich, Phys.Lett.A 233 (1997) 30. 

[12] S. Meljanac, M. Milekovic and A. Samsarov, Phys. Lett. B 573 (2003) 202 
[arXiv:hep-th/0307126]. 



13 



[13] S. Meljanac, M. Milekovic, A. Samsarov and M. Stojic, Mod. Phys. Lett. B 18 
(2004) 603 [arXiv:hep-th/0405132]. 

[14] S. Meljanac and A. Samsarov, Phys. Lett. B 600 (2004) 179 [arXiv:hep- 
th/0408241]. 

[15] T. Guhr and H. Kohler, arXiv:math-ph/0408033. 

[16] A. Jevicki and B. Sakita, Nucl. Phys. B 165 (1980) 511. 

[17] I. Andric and V. Bardek, J. Phys. A 24 (1991) 353. 

[18] I. Andric, V. Bardek and L. Jonke, Phys. Lett. B 357 (1995) 374 [arXiv:hep- 
th/9411136]. 

[19] I. Andric and D. Jurman, Fortsch. Phys. 50 (2002) 675. 

[20] A. P. Polychronakos, Phys. Rev. Lett. 74 (1995) 5153 [arXiv:hep-th/9411054]. 
[21] M.L.Mehta, Random Matrices, 2 nd edition, Academic Press, N.Y. (1991). 



14 



